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SUMMARY 


The  downwash  angle  and  its  rate  of  change  with  foil  angle  are  derived 
for  a hydrofoil  moving  under  a free  surface  disturbed  by  sea  waves.  By  means 
of  a linearized  analysis,  the  downwash  and  its  rate  of  change  are  found  to 
be  the  sum  of  the  values  determined  for  motion  under  an  undisturbed  free 
surface  and  a sinusoidal  time -dependent  component  due  to  the  sea  waves.  The 
phase  angle  and  frequency  relationships  for  the  time -dependent  components 
are  found  for  both  following-sea  and  head-sea  conditions. 
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INTRODUCTION 

In  order  to  study  the  stability  of  a tandem  hydrofoil  system  in 
waves,  it  is  necessary  to  determine  the  forces  acting  at  the  rear  foil 
and  their  derivatives.  The  rear  foil  will  be  in  a flow  field  composed 
of  the  wave  system  through  which  the  configuration  moves  and  of  the  wake 
of  the  forward  foil.  Since  the  force  derivatives  for  the  rear  foil  are 
affected  by  the  downwash  through  inclusion  of  terms  involving  de/da  (rate 
of  change  of  downwash  angle  with  respect  to  angle  of  attack  of  the  for- 
ward foil)  in  the  equations  of  motion  (Reference  1),  the  necessity  of  de- 
termining the  downwash  pattern  is  evident. 

Ir.  the  present  report,  an  analysis  is  made  of  the  effect  of  ocean 
waves  on  the  lift  force  of  a single  hydrofoil  restrained  in  heave  and 
pitch,  based  on  the  classical  wave  theory  which  assumes  a sinusoidal 
wave  form.  The  resultant  wake  disturbance  generated  by  the  hydrofoil 
(time-dependent  at  any  point  located  at  a constant  distance  aft  of  the 
foil),  the  downwash  angle  £,  mid  de/da  are  then  determined.  The  result- 
ant water  surface  shape  at  any  distance  aft  of  the  foil  is  also  found 
in  this  analysis,  giving  a representation  of  the  wave  pattern  in  which 
the  second  foil  of  a tandem  system  will  be  operating. 

The  study  was  carried  out  at  the  Experimental  Towing  Tank,  Stevens 
Institute  of  Technology,  under  Office  of  Naval  Research  Contract  No. 

Nonr  263-01. 
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A = aspect  ratio 

a - surface  wave  amplitude 

C,  = quasi-steadv  lift  coefficient 

C,  = total  lift  coefficient  in  waves 
' lt 

- quasi-steady  two-dimensional  lift  coefficient 
c = wave  velocity 

c'  = mean  hydrofoil  chord 

g = acceleration  due  to  gravity 

H = vertical  displacement  of  hydrofoil  from  equilibrium,  positive  down- 
wards 

h = submergence  of  foil  below  smooth  water  surface 

s = distance  aft  of  1/4-chord  point  of  hydrofoil 

t = time 

u0  = longitudinal  component  of  orbital  velocity 

V = velocity  of  foil 

wQ  = vertical  component  of  orbital  velocity 

wQ  - v>0/V,  change  in  angle  of  attack  due  to  orbital  velocity 

a = hydrofoil  angle  of  attack 

7 = 7'  = circular  frequency  of  encounter  in  following  sea  and  head  sea,  re- 
spectively 

e = downwash  angle 

£s  = disturbance  wave  produced  in  smooth  water 
= disturbance  wave  produced  in  waves 
77  = sea  wave  surface  displacement 

7?r  = resultant  wave  surface  displacement  aft  of  hydrofoil 

\ = wave  length  of  sea  waves 

<£,< i/,8, o,9  = phase  angles 

Subscript  e denotes  equilibrium  condition 
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BASIC  ASSUMPTIONS  UNDERLYING  THE  THEORY 

This  study  is  limited  to  the  motion  of  a two-dimensional  hydrofoil 
restrained  in  heave  and  pitch  and  moving  at  constant  speed  under  a free 
surface  disturbed  by  sea  waves.  A linearized  treatment  is  utilized  in 
considering  the  total  lift  to  be  the  sum  of  the  quasi -steady  lift  due  to 
motion  at  a fixed  geometric  angle  of  attack,  the  lift  variation  due  to 
the  orbital  velocities  of  the  waves,  and  the  lift  variation  due  to  change 
in  submergence,  within  the  limits  of  the  region  of  stall . The  forces  are 
assumed  to  respond  instantaneously  to  the  wave  influence.  No  consideration 
is  given  to  the  effect  on  the  foil  lift  of  the  vortices  cast  off  ir.  the 
wake  due  to  the  nonstationary  flow.  This  assumption  is  based  on  the  fact 
that  the  frequency  is  low  and  is  not  comparable  in  magnitude  to  flutter 
frequencies.  Thus,  for  stability  analysis  purposes,  the  unsteady  flow  ef- 
fects can  be  neglected. 
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PROPERTIES  OF  THE  WAV2  MOTION 


The  wave  shape  for  a simple  harmonic  progressive  wave  moving  from 
left  to  right,  and  defined  with  respect  to  axes  fixed  in  the  fluid,  is 
given  by 

2v 

V = a sin  x (x  - ct ) (1) 

The  origin  of  coordinates  A'  - 0,  1'  = 0 is  taken  at  the  undisturbed  water 
surface,  and  t = 0 is  taken  at  the  instant  the  nodal  point  passes  the 
origin.  The  distance  A is  positive  in  the  right-hand  direction  and  Y is 
positive  upwards. 

In  a progressive  harmonic  wave,  the  fluid  particles  move  in  cir- 
cular orbits  and  the  horizontal  and  vertical  components  of  orbital  ve- 
locity and  orbital  acceleration  at  a depth  h are 


uo  = ° (nr)  sin-X-(*  - ct) 

"o  = -«  (nr)  Cos-XT*  “ ct)  , 

“o  = -a  (nr)  «'2w*/A  cosx-fA  - ct)  , 

io  m~a  (rr)  *~ivk/K  sin^-fA  - Ct)  . 


cos "(X  - ct ) , 


9tt 

sin-T—fA  - ct) 


The  absolute  magnitude  of  the  orbital  velocity  is  then 


- (-¥-) 


For  a wave  propagating  from  right  to  left,  the  above  equations  for  wave 
profile,  orbital  velocities,  and  orbital  accelerations  are  valid,  pr<  - 
vided  the  sign  of  the  wave  profile  velocity,  r,  is  changed. 
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TOTAL  LIFT  VARIATION  IN  WAVES 

The  total  lift  developed  by  the  hydrofoil  is  equal  to  the  sum  of  the 
quasi-steady  lift  of  the  foil  moving  under  an  undisturbed  free  surface,  the 
change  in  lift  due  to  orbital  motion,  and  the  change  in  lift  due  to  the 
change  i rgence  c£  the  foil. 

The  quasi-steady  lift  coefficient  for  a two-dimensional  hydrofoil  is 
given  in  Reference  2 as 

cJLs=  2 mg(Px  - a.gP2)  , (5) 


where  ag  is  the  geometric  setting  relative  to  the  smooth  water  surface 
(two-dimensional  angle  of  attack),  and  P,  and  P,  are  dimensionless  ex- 
pressions dependent  on  depth  and  Froude  number  V2/gc'  , which  express  the 
effect  of  the  proximity  of  che  free;  surface.  For  the  small  angles  usu- 
ally experienced  in  the  motion,  the  two-dimensional  lift  coefficient  may 
be  approximated  by 

_ C£s  = 27TCLgPl  . (6) 

For  the  purpose  of  applying  the  results  of  this  investigation  to 
stability  studies,  a more  reliable  effect  can  be  expected  by  utilising 
the  expressions  for  lift  coefficient  and  its  derivatives  that  are  given 
by  three -dimensional  considerations.  Thus,  the  resulting  wake  disturbance 
will  then  be  connected  with  the  forces  produced  by  an  actual  three-dimen- 
sional hydrofoil.  For  an  actual  three-dimensional  foil,  it  has  been  de- 
termined in- Reference  3,  by  means  of  the  momentum  theory  for  airfoils, 
that 


ac, 



da 

and  hence 


2 irP1- 


*rrA  J.  Q U / 

• w-/  v 


774  + (8h/c’)  + 4ttP1 


(7) 


if 


",  = 2vaP. 

Ls  1 77 A + 


77/1  + 8 h/c' 


v A + (&h/c‘  ) + 4^1 


/ n \ 

VO/ 


where  a in  this  case  is  the  actual  angle  of  attack  of  the  foil  relative  to 
the  zero  lift  angle. 
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The  form  of  CL  given  by  equation  (8)  may  now  be  used  to  determine  the 
derivative  of  lift  coefficient  with  respect  to  depth,  which  is  found  to  be 


ifb 


2 tra 


■nA  + 8 /*  / c ' 


3P, 


, Ci2-nPJc‘  ) - 477  ^-  (vA  * Bh/c‘  ) 

+ p. — 

nA  + f8Vc'  ) + 4^!  [,7(4  + (Qh/c’  ) + 477P  ] 2 


(9) 


Since,  in  the  linearized  analysis,  the  coefficients  are  to  be  evaluated  at 
the  equilibrium  position  and  since  the  definition  of  H is  deviation  from 
equilibrium  depth,  the  relation 


(10) 


readily  follows.  The  above  derivations  may  also  be  found  in  Reference  1. 

The  major  effect  of  the  orbital  motion  on  the  lift  of  a hydrofoil  is  to 
change  the  angle  of  attack.  The  angular  change  in  linearized  form  is  expressed 
as 


A 


(ID 


where  V is  the  forward  speed  of  the  hydrofoil.  By  ignoring  the  lift  variation 
due  to  orbital  accelerations  and  the  change  in  longitudinal  speed  (orbital 
velocity  uQ),  the  lift  variation  due  to  orbital  motion  may  be  approximated  by 


(12) 


as  shown  in  Reference 


The  change  in  submergence  occurs  because  of  the  local  elevation  or  de- 
pression of  the  fluid  surface  due  to  the  wave  motion.  The  change  in  lift  co- 
efficient due  to  change  in  submergence  of  the  foil  is  then  -'.pressed  a.s 


/ac,  \ 

hsV 


(13) 


The  total  lift  coefficient  for  the  hydrofoil  in  waves  is  now  given  by 

(14) 


s - c-.  *&)*•♦&)  ^ • 
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WAKE  DISTURBANCE  PRODUCED  BY  THE  HYDROFOIL 


In  studying  the  flow  about  a body  moving  with  velocity  V with  respect 
fo  the  undisturbed  fluid,  it  is  more  convenient  to  refer  to  axes  fixed  in 
the  body  rather  than  in  the  fluid.  The  required  transformation  is 


* = X - Vt 

y - Y 

t = t 


Y (15) 


where  the  origin  x = 0 is  taken  at  the  quarter -chord  point  of  the  hydrofoil 
and  x is  positive  to  the  right  (in  the  direction  of  motion).  The  wave  pro- 
file and  orbital  vertical  velocity  are  then  given  as 


V * a sin + (V  - c)t) 

A 

*°  * cos  Tf*  *"  (v  ~ c)* i 

for  a following  sea  and  as 


Y (16) 


V m a sin  -£”[*  + (V  * c)t) 
A 


V 


o 


cos  + (V  + c)t] 


Y (17) 


for  a head  sea.  The  frequency  of  encounter  for  following  seas  and  head  seas 
is  given  by  (V  - c)Aand  (V  + c )A,  respectively. 

For  • hydrofoil  aoving  at  speed  V and  situated  at  depth  k under  an 
undisturbed  free  surface,  the  resultant  surface  wave  disturbance  (beyond 
the  first  1/4  wave  length  downstrean)  is  given  in  Reference  5 as 


--C, 


** 

sin — 7 
Vt 


(18) 


where  t is  the  distance  neasured  aft  of  the  hydrofoil  (s  m-x).  Since  the 
wave  disturbance  is  proportional  to  the  foil  lift  coefficient,  it  has  been 
assuaed  that  the  sane  effect  is  true  for  the  foil  in  waves  where  the  lift 
has  a tine -dependent  sinusoidal  variation.  Thus,  the  wake  disturbance  pro- 
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duced  by  the  hydrofoil  operating  in  waves  will  also  have  a time-dependent 
sinusoidal  variation  and  is  given  by 


^ " ^Cltc' e't>,/v2  sin-^ 


= - C.  ♦ 


. + V 


, n\\ 

+ Vah^T 


*e-«‘/r2sin-^- 

V2 


Substituting  the  expressions  for  »0  and  17  from  equations  (11),  (16),  and 
(17)  gives  the  total  lift  coefficient  in  a following  sea: 


lt  " c*.  + °v(ir)  + fe)  e'4”k/A(~*r)  sin[Trv  • c;‘ " +_ 


where 


<t>  = tan 


-bV,nw&) 


There  is  no  dependence  on  the  coordinate  x since  the  lift  ac*"  at  the 
quarter -chord  point  of  the  foil,  where  x 3 0.  Equation  (20)  may  be  simpli- 
fied by  letting 


b .praws 

V\Bh/  W/  \ 9a/ 


277  i 

7- T"-"  J 

where  y is  now  a circular  frequency  of  encounter,  resulting  in  the  form 

C,  - C,  + a B sin (yt  - <t>)  ( 

Lr  s 

Similarly,  for  the  case  of  a head  sea, 

C.  m C + ag  sinfy' t + <t>)  , ( 

*■  r 


Security  information 
CONFIDENTIAL 


[ 


R-497 

-10- 


COHf IDENTi AL 
Security  Information 


where 


y'  - ~(V  ♦ c)  (25) 

is  the  circular  frequency  of  encounter  for  head  seas. 

The  wake  disturbance  produced  by  the  hydrofoil  operating  in  waves  is 

then 

**  -fc^  + °B  sin^t  “ ^)Jc'e'*k/r2  sin-py  (26) 

for  a following  sea,  anu 

C.  - -[ Ct<  + afi  sinfy't  + d>)Jc' e-«fc/rJ  sin-^-  (27) 

for  a head  sea. 

With  the  mathematical  development  of  the  disturbance  wave  completed, 
it  is  now  possible  to  determine  the  water  surface  shape  at  any  distance  aft 
of  the  foil,  which  will  then  represent  the  wave  pattern  in  which  the  second 
foil  of  a tandem  system  will  be  operating.  At  a distance  JL  aft  of  the  quar- 
ter-chord point  of  the  foil  (x  * -X)  , the  original  sea  wave  pattern  will 
be  represented,  st  time  t,  by 


i 

H 

I 

1 

■ 

■ 

■ 

■ 

■ 


7j  ■ a sin(>t  - i p) 
for  a following  sea,  and  by 
17  ■ a sin (y't  - \p) 
for  a head  sea,  where 

2-nl 


(28) 


(29) 


4> 


(30) 


However,  interest  is  centered  on  the  resultant  wave  at  the  time  ( + J l/V , 
since  at  this  time  the  disturbance  created  Ly  the  foil  at  time  t is  felt 
U'jin  the  rear  (*  * -X)  , because  of  the  time  lag  effect.  Thus,  st  time 
* ♦ i'/V,  the  resultant  wave  surface  st  x * - X in  s following  sea  is 
represented  by 
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■ 


. / y L 

Vr  = a sin  Vt  + \b 

\ V 


+ r 


(31) 


where  is  the  wave  amplitude  at  s = II,  due  to  a disturbance  produced  at 
time  t.  In  this  derivation,  it  is  implicitly  assumed  that  the  existing  sea 
wave  pattern  and  the  disturbance  wave  are  superposable  and  that  no  inter- 
ference exists  between  them.  The  resultant  wave  pattern  at  * - — JL  in  a 
following  sea  may  now  be  expressed  as 


Vr  - C sinfyt  + 8 ) - C c' e'*k/v  sin 

**  c 


(32) 


where 


= <rWl  +B2c'  2e'2*h/v  sin2  - ZBt  • e *k/v  sin  -py  cos^^  -</»  + 


8 = tan 


-1 


inf—  - i ji\  + ftc‘ e~*k/v  sin -^4  si 


nd> 


cos 


(33) 


if—  - *p\  - Be’  e'**/k  sin-^y  cos<£ 

V I v 2 

Similarly,  the  resultant  wave  shape  at  x =>-JZ,  ir.  s head  sea  is  expressed  as 


2 op 

vr  = o sinfy't  » crj  - c.  c'  e~tk/v  sin — =- 

‘s  v 


(34) 


where 


<^/l  + B2c'  sin2  -^y  -26c'  e-**/K2  sin  -^y  cosf2^  -0 -<^ 


(35) 


O'  » tan 


sinf— — - - Be'  e~*h/y  sin-^4-  sind> 

\ V / ^ 


cos 


\ r 1 


^ 2 C 

^ ) - Be' e',*/K  sin — i- cosd> 
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DETERMINATION  OF  DOWNWASH  ANGLE  £ AND  OF  de/da 


The  angle  of  attack  at  the  rear  foil  of  a tandem  hydrofoil  system  is 
dependent  upon  the  mean  downwash  angle,  £.  The  sense  of  £ is  positive  when 
it  tends  to  reduce  the  angle  of  attack  at  the  rear  foil.  It  is  determined 
from  the  slope  of  the  sinusoidal  disturbance  wave  generated  by  the  bound 
vortex  of  the  forward  foil.  The  downwash  angle  £ is  the  flow  angle  at  the 
rear  foil  due  to  the  disturbance  wave  created  by  the  forward  foil.  Thus, 
it  is  not  to  be  considered  as  being  due  to  the  influence  of  the  orbital 
velocity  of  the  sea  wave  profile  above  the  rear  foil;  it  is  found  from  the 
expression  for  and  not  vr • 


For  a hydrofoil  in  waves,  the  disturbance  wave  is  represented  by 
equations  (26)  and  (27)  for  following  seas  and  head  seas,  respectively. 
The  wave  slope  at  the  surface  is  then  given  as  d^/ds;  the  angle  of  down- 
wash  at  a distance  « behind  the  foil  and  at  a depth  h'  below  the  water 
aurface  is  expressed  as 


...„Z 

£ ■ g”f  * ^ 

ds 


[c 

L ‘i 


♦ =B  sinfyt 


COS 


g® 


(36) 


for  a following  sea,  and  by 


♦ aB  sin (y‘  t + <t>) 


+ CQS 


(37) 


for  a head  sea. 


By  using  the  preceding  expressions  for  £,  the  rate  of  change  of  down- 
wash,  de/da,  can  be  easily  found.  The  terms  in  the  expressions  for  £ that 

are  functions  of  a are  C. , B,  and  <t>:  the  derivative  9C,  /9a  is  given  by 

T*  . .* 

equation  (7)  and  the  derivatives  9g/9a  and  90/9a  are  given  below  as 


9B_  _ 5 (*»,Y 

9a  26  9a  \ 9 h j 


4nJa 


vA  * 8 h/c' 


9 P, 


(32TTPJ/C'  ) -4v—±-(v a +8 h/c' ) 
oh 


B YnA  * (Bh/c1  ) *ivPl  9h 


♦P, 


[vA  + (8 h/c'  ) + 4^]* 


(38) 


I 
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tan  <t>  3 ( Ls 


-tan  4> 


1 + tan ii>  l^Lsj  (1  + 


tanJ0)a 


The  value  of  u’c/u'u.  for  a hydrofoil  in  a following  sea  is  then 


Llf-'NA®  * eVJ  ^ 


where 


^ ’•  tan-1 


Similarly,  for  the  hydrofoil  in  a head  sea, 


co.-^l 


The  change  in  angle  of  attack  at  the  rear  foil  of  a restrained  tandem 
hydrofoil  system  in  waves  is  given  by 

Aa  » »c  - c , (43) 

where  »0  and  e are  evaluated  at  the  location  of  the  foil  and  are  functions  of 
the  wave  properties. 
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CONCLUSIONS 

The  equations  developed  in  the  present  report  give  expressions  for 
tne  downwavh  agle  and  its  rate  of  change  with  foil  angle,  for  a hydrofoil 
moving  under  a free  surface  disturbed  by  sea  waves.  The  total  wave  disturb- 
ance is  assuned  to  be  a linear  superposition  of  the  effects  of  the  quasi  - 
steady  lift  and  the  time -dependent  sinusoidal  lift  variation  due  to  orbital 
motion  and  cnange  of  submergence.  The  resultant  downwash  angle,  t,  and  its 
rate  of  change,  de/da,  are  then  the  sum  of  the  contribution  found  for  the 
case  of  motion  under  an  undisturbed  free  surface  (Reference  1)  and  a sinu- 
soidal time-dependent  component  due  to  the  sea  wave  influence.  The  phase 
angles  for  the  time -dependent  components  of  e and  de/da  (relative  to  the 
phase  of  the  disturbing  sea  waves)  arc  negatives  of  each  other  for  the  fol- 
lowing-sea and  head-sea  conditions.  The  frequencies  for  the  time -dependent 
terms  in  e and  de/da  are,  however,  fundamentally  different  for  the  following 
sea  and  head  sea. 
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